Scalar Wave Propagation in Random, Amplifying Media: Influence of Localization 
Effects on Length and Time Scales and Threshold Behavior 
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We present a detailed discussion of scalar wave propagation and light intensity transport in three 
dimensional random dielectric media with optical gain. The intrinsic length and time scales of such 
amplifying systems are studied and comprehensively discussed as well as the threshold characteristics 
of single- and two-particle propagators. Our semi-analytical theory is based on a self-consistent 
Cooperon resummation, representing the repeated self-interference, and incorporates as well optical 
gain and absorption, modeled in a semi-analytical way by a finite imaginary part of the dielectric 
function. Energy conservation in terms of a generalized Ward identity is taken into account. 

PACS numbers: 



I. INTRODUCTION 

Built on a wide story of success and plenty of achieve- 
ments in science, light propagation and intensity trans- 
port in multiple-light-scattering random media further 
provides a lot of remarkable and fascinating features in 
theory^— and experiment^. Special ongoing interest is 
focused on optical gain material regarding disordered as 
well as periodic structures^—. 

In theoretical physics a profound understanding of this 
subject requires the discussion of strong or Anderson 
localization 12 of light, a microscopic transport theory in 
the diffusive limit usually based on conservation laws, 
incorporation of gain or absorption into generalized con- 
servation laws, as well as observing the occurrence of a 
threshold behavior of the light intensity. 



Anderson localization has been show] 



.13.14 



to originate 



from repeated self-interference of diffusive modes. In 
experiment*^— the authors have tested this. Since in- 
tensity diffusion is solely based on energy conservation, 
a rigorous and consistent framework is needed to de- 
scribe the interplay between coherent amplification and 
localization. Especially since coherent amplification is 
expected to enhance transmission whereas localization 
tends to stop light transport altogether. This interest- 
ing subject is also discussed in the context of random 
lasingiiS, where it has been shown by measuring the pho- 
ton statistics 2 ^, that the laser emission is due to coher- 
ent feedback and occurs from spatially confined spots in 
the sample. A theoretical attempt^ to independently 
explain such phenomena has been proposed based on 
scatterers statistically forming ring resonators within the 
sample, which is able to provide a feedback mechanism. 

Despite the issue of transport and localization, also the 
onset of lasing behavior, the crossing of the so-called laser 
threshold, has triggered still ongoing publications^^ 2 - - — 
for the last decades. On the theoretical side, the difficulty 
lies in the use of static gain below threshold, which is per- 
fectly valid unless the threshold is reached. Beyond the 
laser threshold, the laser dynamics have to be explicitely 



taken into account. Such a threshold behavior has to be 
carefully respected and incorporated in a consistent way 
into any theory of transport and localization. 

In the present paper we clarify this subject by pre- 
senting a semi-analytical general theory of light propaga- 
tion, localization effects in the presence of optical gain, 
and consistently discuss the occurrence and influence of 
a growth threshold on the single particle propagator as 
well as on the energy density correlation function, related 
to the intensity in the sample. In this way we set the last 
stage within the range of linear response theory, which 
then will serve as a basis to involve the actual lasing dy- 
namics. This dynamical behavior is, however, not sub- 
ject of this article and will be presented in forthcoming 
publication. 



II. MODEL AND THEORY 

A. Basic Setup 

Systems of significant experimental relevance^— con- 
sist of (almost) spherical scatterers embedded into a 
background medium forming some emulsion. For a theo- 
retical description we therefore consider identical spher- 
ical scatterers located at random positions. The scatter- 
ers as well as the background medium are respectively 
assumed to be homogeneous and hence will be described 
by dielectric constants e s and e^, respectively. Within a 
scmiclassical context linear absorption and optical gain 
will be represented by a finite positive or negative part of 
the dielectric function, so in general Ime s =/= ^= line?, is 
assumed. Throughout the paper we neglect polarization 
effects and therefore consider the scalar wave equation 
which has been Fourier transformed from time t to fre- 
quency u> and reads 

^ e(r )* u (?) + V 2 V u (r) = -iu>^-j u (r) , (1) 

where c denotes the vacuum speed of light and j u (r) the 
current. The dielectric constant e(r ) = + AeV(r), 
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where the dielectric contrast has been defined as Ae = 
e s — tb, describes the arrangement of scatterers through 
the function V(r) = (r— R), with (r) a local- 

ized shape function at random locations R. The inten- 
sity is then related to the ficld-ficld-corrclation function 
(^(f, t )^*(r', t ' )) here angular brackets (. . .) denote en- 
semble or disorder average. To calculate the field-field- 
correlation the Green's function formalism is best suited, 
the (single-particle) Green's function is related to the 
(scalar) electrical field by 



W, * ) = Jd 3 r ' Jdt ' G(f, f';t, t> )j(f' , t> [ 



(2) 



The Fourier transform of the retarded, disorder averaged 
single-particle Green's function of Eq. (fTJ) reads, 



/~fU 

* ~ e b {u/cY - \qf 



^q 



(3) 



where the retarded self-energy arises from scattering 
off the random "potential" — (uj/c) 2 (e s — €b)V(r ). Using 
the Green's function the mode density N(uj) may be ex- 
pressed as N(u>) = — (w/7r)ImGo , with the abbreviation 
used throughout this publication G% = J d 3 q/(2ir) 3 GZ. 

In order to study the transport of the above intro- 
duced field-field-correlation we consider the so-called 4- 
point correlation function, defined in terms of the non- 
averaged Green's functions G, G* in momentum and fre- 
quency space as (Q,fi) = (^ + + |-^>. Here we 

have introduced the usual 3 ^ center-of-mass (q, oj) and rel- 
ative (Q, fl) frequencies and momenta: The variables fi, 
Q are associated with the time and position dependence 
of the averaged energy density, with Q = Q/\Q\, while 
uj± = uj ± Q/2 and q± = q± Q/2 etc. are the frequencies 
and momenta of in- and out-going waves, respectively. 

The intensity correlation, or disorder averaged 
particle-hole Green's function, $~,(<3,fi) obeys the so- 
called Bethe-Salpeter equation 



^qq> 



Gf + ( W+ )G^_( W _) 



6{q- ql) 



dV 



q q 



(4) 



By utilizing the known averaged single particle Green's 
function, c.f. Eq. <j3j> , on the left-hand side of Eq. (|4]) 
the Bethe-Salpeter equation may be rewritten as kinetic 
equation, see, e.g., Refi 3 ^, 



(jjQ, 



Rceh 



(q ■ Q) 



-iImG£ 



6(q - qi) + 



qq 



J (27T 



'qq <f <f 



.(5) 



In order to analyze the correlation function's long-time 
(il —> 0) and long-distance (\Q\ —> 0) behavior, terms of 
0(fl 2 , Q 3 , £IQ) have been neglected here and throughout 



FIG. 1: Ladder approximation of the total particle-hole ver- 
tex. The diagrams on the left-hand side form a geometrical 
series and may therefore easily be summed up analytically. 



this paper. Eq. ([5]) contains both, the total quadratic mo- 
mentum relaxation rate 1/t 2 = c 2 lm(ebUJ 2 / c 2 — £") (due 
to absorption/gain in the background medium as well as 
impurity scattering) and the irreducible two-particle ver- 
tex function 7^.,(Q, Q). 

To solve this equation, the technique of expansion into 
moments is used. The technical details of this expansion 
are discussed in the following subsection B. The Reader 
not interested in such details may skip this and readily 
proceed to subsection C. 

Furthermore it is to be noted that the energy con- 
servation is implemented into the solution of the Bethc- 
Salpctcr equation in a field theoretical sense by a Ward 
identity (WI) which has been derived for the photonic 
case in Ref. l36l . and which for scalar waves takes the ex- 
act form 



(d 3 q ' 
" J(2tt) 

= /«(fi) 



q+ q- 



q+ Q- 



ReS 



^If-ffftSi) (6) 



The right-hand side of Eq. ([6]) represents reactive effects 
(real parts), originating from the explicit w 2 -dependencc 
of the photonic random "potential" . In conserving me- 
dia (Imet, = Ime s = 0) these terms renormalize the 
energy transport velocity v E relative to the average 
phase velocity c p without destroying the diffusive long- 
time behaviori 27 ! 36 In presence of loss or gain, however, 
these effects are enhanced via the prefactor = 
(cj$lRcAe + iw 2 ImAe)/(w 2 ReAe-Hw^ImAe), which now 
does not vanish in the limit Q — > 0. 



Expansion of Two-particle Green's Function 
into Moments 



In order to extract a diffusion pole structure out of 
the Bethe-Salpeter equation Eq. (H}, the correlator or 
cquivalently the q ' integrated correlator 



(7) 



has to be decoupled from the momentum dependent pref- 
actors with the help of some approximation scheme. In 
this subsection we discuss this procedure in analogy to 
the argumentation for electronic correlations presented 
in reference^. Such an approximation must obey the 
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results of the so-called ladder approximation as well as 
it must incorporate the set of physical relevant variables 
involved in observed phenomena. In a first step we use 
the bare first two moments of the correlation function ^ 



defined according to 



<f> jp (Q,Q) 



( d 3 q f dV 
J (2tt) 3 J (2tt) 3 

r d 3 9 r dV 

J (2tt) 3 J (2tt) 3 



(8) 
(9) 



respectively. The second step is to recognize that these 
bare moments are related to physical quantities, the 
energy density correlation P£(Q,SY) and the current- 
density-correlation J£(Q, £1), by dimensional prefactors: 



to 

UJV 



UJV E 



P£(Q,n) (io) 
j-(Q,n). (ii) 



The projection of the correlator Eq. (J7)), onto the 



bare 



moments $ PP (<f, fl) as defined in Eq. 



■j P {q, fl), shown in Eq. (|9|), is therefore given by 
fdV^ A(C 
J( 



d 3 g' 

(2tt) 3 



<MQ,^) 



g(g)(g • g) 



and 



(12) 



*ip(Q,fi). 



where the projection coefficients A(<f ) and i?(cf ) are to be 
determined in the following. For obvious reasons in this 
expansion the bare moments may be substituted by their 
physical counterparts energy density P£ in Eq. (fTU| and 
current density J" (Q, H) from Eq. ([TT|) . The expansion 
coefficients A{q) and B{q) in Eq. (fT2j) behave uncriti- 
cally under localization, so they can be determined using 
the simple ladder approximation, where all expressions 
are known exactly. The ladder approximation of the two- 
particle vertex function is explicitely illustrated in Fig. [T] 
In the following we use this approximation and demon- 
strate how to obtain the expansion coefficients from it. 
In ladder approximation the zeroth bare moment is given 
by: 



r d 3 g 

J (2tt> 
7o 

the superscript i 
proximation and in 



G 9 - + (Q,fi)G*-_(Q,0) 



(13) 



refers to 
the last 
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the 
step 



ladder ap- 
the product 



has been expanded 

up to linear order in q furthermore the renormalizcd 
vertex 70 is given by 

(Re 70 G + Re E) w 2 Ime b 



7o = To+iU^)- 



Im Go 



ImGo 



(14) 



where 70 is the bare vertex and /^(fJ) arising from the 
Ward identity has been defined in Eq. ©. Within the 
simple ladder approximation the bare moment $^ (Q, fi) 
defined in Eq. §§§ is thus given by 



*&(Q,n)= 



dV 



' d 3 q _ - 

7^33(9 ■ Q) G q+ G }-\ j2^ G q'+ G *q' T L -(15) 



Following the above strategy and expanding the product 
G^i Gt, under the second integral up to first order in q 

one obtains the expression 

«*(&n) = lr J^or • q)g 9 - + g;-_ (ie) 

By now employing the same expansion to the remaining 
product of the Green's function one eventually finds 



T L f d 3 <z . l AGjg-QjQ 
^J(2^ (g ' Q) - — ' (17) 



70AG0 



where the abbreviation AG = G — G* has been intro- 
duced and will be used throughout this paper. 

In the next step of determining the expansion coeffi- 
cients A(q) and B(q ) defined in Eq. ([12")) we go back 
to the field-field correlation function . Within the 
uncritical ladder approximation the two particle Green's 
function is given by 

f dV 



J(2tt 



Employing again the momentum expansion of the single 
particle Green's function the above equation, Eq. (fT8l 
can be simplified to yield 



<I>„ 



AGj- i AG}(g-Q)Q 
7o 2 AGo ii+ 2 7o 2 AG L ' 



(19) 



Finally we are in the position to start putting things to- 
gether. By using the above given momentum expansion, 
Eq. (fT9|) . together with the expressions given in Eq. (jTTJ) 
and in Eq. (|13[) in conjunction with the proposed projec- 
tion, or expansion into moments, Eq. (|12l) . the following 
relation is eventually obtained 



AG,- l AG}(q-Q)Q 
7 2 AGo L+ 2 7o 2 AGo 1 



(20) 



A(q] 



1 



£?(?)(?• Q) 



J^£(gO(?'-Q) 2 ^ 



70 J (2tt) 3 



2 70AG0 



By comparison of coefficients in the above relation, Eq. 
(|20| . the demanded coefficients A{q) and 5(g) of the ex- 
pansion into moments, Eq. (|12p . can now be determined 
to be 



A(q) 



AG f 



5(<f) 



AG}. 



(21) 
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Employing those expressions for the expansion coeffi- 
cients, one may eventually express the two-particle cor- 
relator •f'gVf in the following way 



I 



(2tt)3 f? ' 



-TO,fi) 



(22) 



AG, 



9' 



AG|(g ■ g) 



(^)/&AG|,(r-o) : 



<(Q,n). 



The above expression, Eq. (|22|) . represents the complete 
expansion of the intensity correlator into its moments. 
This will be used in the next subsection to decouple and 
therefore solve the Bethe-Salpeter equation. 



C. General Solution of the Bethe-Salpeter 
Equation 

The disorder averaged intensity correlation, the two- 
particle Green's function, obeys the Bcthc-Salpctcr equa- 
tion, see Eq. 



Q 1 q+ q- 



fdV 



Iqq" ^q"q> 



(23) 



as already discussed the Bethe-Salpeter equation may be 
rewritten into the kinetic equation given in Eq. (JS|) 



^q 



wft2Ree - Q ( q ■ Q ) + AS - u/Ae 
fdV 

J(2^ AGf7 ^^'- (24) 



To find the solution of Eq. (|24j) . in a first step one sums 
in Eq. (|24j) over momenta <f, incorporates the general- 
ized Ward identity as given in Eq. and subsequently 
expands the obtained result for small internal momenta 
Q and internal frequencies f2. It is also essential to em- 
ploy the decoupling shown in Eq. (|2"2"|) . Eventually after 
some algebraic manipulations the generalized continuity 
equation for the energy density is found to be 



np^ + q j e " = 



Alii uj N(oj) 



| »[gi 0) +AM] 



(25) 



which represents energy conservation in the presence of 
optical gain and/or absorption. 

Within the standard solution procedure the next step 
is to obtain a linearly independent equation which also 
relates the energy density P£ and the current density J£ . 
This is realized in a similar way to above, foregoing one 
first multiplies the kinetic equation, Eq. (|24|) . by the pro- 
jector q ■ Q and then follows the above outlined recipe 



to eventually obtain the wanted second relation, this is 
the so-called current relaxation equation 



Ree;, 



J£ + ^ QP£ = , (26) 



relating as demanded energy density P£ and energy den- 
sity current J£ and furthermore introduces the so-called 
memory function M(Q) according to 

M(n) = ^I^W-Ql^'m^m'-Ql (27) 



d^q 
(2^F 



[?-Q] 2 (A^) ! 



where 7—, = ypgi(Q, ^) is the total irreducible two- 
particle vertex, which will be discussed in more detail 
in the following subsection. 

So far, two independent equations, Eq. (f2"5]l and Eq. 
(|26|) . have been obtained, both of them relating the cur- 
rent density and density P£ . Therefore one may now 
eliminate one of the two variables in this linear system of 
equations. One chooses to combine the two equation to 
find an expression for the energy density 



4™7VH/( 5 £ 1) [1 + A(w)]4) 



n + iQ 2 D 



(28) 



exhibiting the expected diffusion pole structure for non- 
conserving media, i.e. in the denominator of Eq. (f2"5)) 
there appears an additional term as compared to the case 
of conserving media. This is the term £~ 2 D, sometimes 
referred to as the mass term, accounting for loss (or gain) 
to the intensity not being due to diffusive relaxation. In 
Eq. |28p also the generalized, Independent diffusion co- 
efficient D(il) has been introduced via the relation 



£>(fi) [I - i nujT 2 Ree b ] = D' ° 



2 D(Q)M{uj). (29) 



Furthermore, Eq. ([28]) also introduces the absorption or 
gain induced growth or absorption scale £ Q of the diffusive 
modes, 



f 



r e A e — 2w 2 Ime6 



2Ree b - A e B e /uj wD(fi) ' 



(30) 



which is to be distinguished from the single- particle or 
amplitude absorption or amplification length. The dif- 
fusion constant without memory effects in Eq. (|29p . 
Z?q°* = Dq + Df, + D s , consists of the bare diffusion 
constant^, 



_ 2v E c p ( d 3 i 
° ~ ttN(uj) J(2tt 



3 [?.QYQmG% 



w\2 



(31) 



and renormalizations from absorption or gain in the back- 
ground medium (D b ) and in the scatterers (D s ), 

D b = {ojt) 2 lme b Do/4 , D s = r £ A e T 2 D /8 , (32) 



Y"(Q,S!> = >< 
kk 



k + E, 



k + k' + 



k_ k'_ 



FIG. 2: The upper panel shows a diagrammatic expansion 
of the irreducible two-particle vertex 7. The lower panel dis- 
plays the disentangled Cooperon with changed momentum 
arguments as discussed in the text below. 



£ a 2 D. Finally the memory kernel M(f2) reads 



M(n) 



(2v E c p ) 2 u e [2?ra;u e 7V(a;) + r e A e - 2u 2 lme b ] 



( d 3 q ( d 3 q' 
) (2tt) 3 J (2tt) 3 



ttujN(lu)D D({1) 
[g-Q]\lmG q \{lmG ql ) 2 W ■ Q] 
^ + (q + q') 2 +Ha 2 



(33) 



Eqs. (|29| -(|33 |) constitute the self-consistency equations 
for the diffusion coefficient D(£l) and the growth/decay 
length £ a in presence of absorption or gain. 



III. RESULTS AND DISCUSSION 



A. Diffusion Constant and its Renormalization 



where Dq is the same as in Eq. (f3"Tj) . with (ImG^) 2 re- 
placed by Re(G| 2 ). In the above Eqs. (|5D ]) -([5^ ]l the 
following short-hand notations have been introduced, 



-4, 
B, 



Im(Ae£'") 
Im(AeG^) ' 
2[u e RcG + RcI] ] 
(ReAe) 2 + (ImAe) 2 
2w 2 (ReAe) 2 ' 



ImAe/ReAe, 



D. Vertex Function and Self-consistency 



From equations Eq. (J27J) and Eq. (|29| it is clear 
that the energy density or two-particle function given 
in Eq. (f28|) still depends on the full two-particle vertex 



Jgiff- Before discussing the vertex function, we want to 
briefly recall some arguments concerning dissipation. As 
a simplified argument to better understand the physical 
content of the presented systems, one might consider a 
damped harmonic oscillator. The damping term clearly 
introduces dissipation as it breaks time reversal symme- 
try. However, the time reversed solution is still damped 
with the very same damping constant. This shows that 
the dissipation rate itself is invariant under time reversal, 
which at first sight might sound surprising. 

Bearing this in mind one may carefully analyze the ver- 
tex 7^, - for the self-consistent calculation of M (£7) 14 ' 35 , 
exploiting time reversal symmetry of propagation in the 
active medium. In the long-time limit (tt — > 0) the 
dominant contributions to 7^,- are the same maximally 
crossed diagrams (Cooperons) as for conserving media, 
which may also be disentangled. In Fig. ((2]) the disen- 
tangling of the Cooperon into the regular diffusion ladder 
is demonstrated. The internal momentum argument of 
the disentangled irreducible vertex function in the sec- 
ond line of Fig. ([5]) is replaced by the new momentum 
Q = k + k! . By the described procedure 7^, - now ac- 
quires the absorption (gain)-induced decay (growth) rate 



In Eq. ([29]) the diffusion coefficient D(fl,u}) has been 
shown to consist of a memory induced part already dis- 
cussed in the last section and a part Dg ot (u;), given by the 
so-called bare diffusion constant D defined in Eq. (f3"Tj) 
and additional contributions solely due to a finite gain 
or absorption. These renormalizations of the diffusion 
constant are 

D b = (wrf Ime b A)/4 , D s = r £ A e r 2 D /8 . (34) 

Since we are interested in amplifying systems with neg- 
ligible absorption we now want to discuss setups with a 
finite optical gain coefficient inside the scatterers only, 
which are themselves embedded in a conserving media 
e.g. air. In the background line?, is identical zero and 
therefore also D b . The described systems are of strong 
experimental interest^^ and are still not completely 
understood^. 

Before starting with various examples we want to point 
out that in dealing with optically amplifying media one 
has to choose parameters carefully, guaranteeing that the 
system remains below its laser threshold. This will be 
discussed in detail in a following subsection. The pre- 
sented numerical results utilize parameter sets which do 
show bclow-threshold behavior within the considered fre- 
quency range. In particular we present results for three 
characteristic parameter sets, the setup is an optically 
neutral background medium like air (e;, = 1.0), spherical 
scatterers (filling fraction v = 30%) with three different 
gain strengths (e SCQt = 10.0 - {0, le - 4, le - 2}/). The 
system with purely real dielectric functions, i.e. conserv- 
ing media, serves as a reference system, and the gain 
is either typical (Ime sca t = — le — 4) or rather large in 
magnitude, here Ime scat = — le — 2. 

In the upper panel of Fig. [31 the real part of the diffu- 
sion coefficient for a non-dissipative system (black line) 
is compared to systems exhibiting gain (colored lines). 
The diffusion constant has been renormalized to its bare 
coefficient, cf. Eq. ([3~Tj) . As already discussed, small 
gain as compared to threshold gain disadvantages local- 
ization, whereas with increasing gain also the diffusion 
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FIG. 3: Normalized real part (upper panel) and imaginary 
part (lower panel) of the diffusion coefficient D(Q, = 0) for 
different values of optical gain as indicated as a function of the 
dimensionless frequency. For the displayed frequency range, 
the gain value ist below threshold. 



Q 
at 




FIG. 4: Real part of diffusion coefficient D(0) (upper panel) 
and the bare diffusion coefficient Do (lower panel) for differ- 
ent values of the optical gain as a function of the dimension- 
less frequency. The diffusion coefficients are shown in natural 
units {2-nroc), ro is the scatterers radius and c is the vacuum 
speed of light. The influence of optical gain on diffusion is 
seen to be small. 



increases. For the discussed gain values this effect is in- 
verted within higher resonances, because they are much 
closer to threshold, where gain narrowing has already 
overcome this suppression. Although the effects of gain 
on transport are rather small, the small but finite gain 
introduces a completely new feature, an imaginary part 
of the diffusion constant at zero internal frequency, i.e. 
an imaginary part to the dc diffusion coefficient or like- 
wise to the dc conductivity. The normalized imaginary 
part of D is displayed in Fig. [3l normalized to the bare 
diffusion coefficient D Q . In the next subsection it will 
be shown how this finite imaginary part may give rise to 
intensity oscillation within the sample. 

The additional contribution D s from Eq. (|34|) to the 
diffusion constant provided by the amplifying scatterers 
is presented in Fig. [SJ This contribution is a direct con- 
sequence of the photonic Ward identity, Eq. ©, and 
establishes the conserved energy density. However, for 
reasonable gain values considered here, this correction is 
seen to be orders of magnitude smaller than the diffusion 
constant, therefore its influence on the transport proper- 
ties is only weak. 

For completeness and later use we also display the bare 
diffusion constant Dq and the real part of the full diffu- 
sion coefficient. The bare diffusion constant as shown in 
Fig. U exhibits strong variations as function of frequency 
but only small variation with increasing optical gain. The 
full diffusion coefficient in Fig. 2] follows closely the be- 
havior of the bare diffusion, as already indicated in Fig. 
[3J For very small and also for large frequencies the dif- 
ference is negligible, visible effects arc found within an 
intermediate frequency range only. 



B. Length and Time Scales 

Within disordered systems there exist different length 
or time scales, related to both single and two-particle 
quantities. Additionally, a geometrical mean distance 
between each two scatterers r m can be defined by ^ = 

, where ro is the scatterer's radius and v the fill- 
ing fraction. For a filling fraction of v = 30% this ratio 

becomes — = 2.41. 

ro 

The most important single particle length is the so- 
called scattering mean free path l s defined in the Green's 
function 

G ^ = \ — (35) 

where the imaginary part of the self-energy introduces 
the decay length l s 

q = -V^^Mi) + irr (36) 



2Im( yV + ilm£(w)) 

The decay length may equivalcntly be represented as a 
life time of the corresponding k-mode. 

In Fig. |6] the scattering mean free path is shown 
as a function of frequency. The strong variation with 
frequency is known to be typical for the low density 
approximation^— as used in this publication. The mi- 
nor dependence on gain is in agreement with the previous 
subsection and mainly established in slightly more pro- 
nounced dips and enhanced peaks. 
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FIG. 5: Normalized extra contribution to the diffusion co- 
efficient D a /Do for different values of the optical gain as a 
function of the dimensionless frequency. 




FIG. 6: Displayed here is the single particle scattering mean 
free path l 3 for amplifying scatterers and conserving scatter- 
ers. The difference between the two is seen to be small and 
mainly manifests itself in narrow dips and peaks. 



Before proceeding, we want to emphasize that in case 
of real dielectric constants the scattering mean free path 
sets the scale determining the loss due to scattering out 
of a given k-mode. Whereas in case of gain media the 
originally k-modc experiences also an amplification. In 
this way a competition is established between scattering 
and gain. Once the optical gain is strong enough to com- 
pensate the scattering loss, i.e. ImE(w) = 0, that fact is 
interpreted as the crossing of the laser threshold ' ,8 ' 10 - 19 . 
This particular case consequently defines the range of va- 
lidity of the presented theory, which cannot describe the 
onset of the laser dynamics. The gain coefficient is there- 
fore to be chosen such that the system remains below its 
threshold gain value. This will be discussed in detail in 
subsection C below. 

Let us now return to the discussion of the intensity 
and the scales related to it. The two-particle Green's 
function as given in Eq. (|28l ) contains two obvious scales 
originating solely from finite values of the gain coefficient. 
These length scales may be defined by 



4 = 

p — 



2/T 



2tt 



(38) 
(39) 



where £ a represents the amplification or absorption 
length of the intensity and £ osc marks the length over 
which the intensity oscillates, where ^ has already been 
defined in Eq. (|30l ). The corresponding time scales may 
then be defined as 



1 



1 



D 

So 



(40) 

= Q z ImD (41) 

7~osc 

The amplification or growth length £ a is displayed in Fig. 



[7] as a function of the external light frequency. The sin- 
gle scatterer Mie resonances are clearly visible as well 
as strong dependence on the gain value. However, even 
for the strongest presented gain, the magnitude of the 
amplification length remains at least an order of magni- 
tude larger than the corresponding scattering mean free 
path, shown in Fig. H>] Additionally, we have plotted the 
oscillation length £ osc in Fig. [7] As compared to the am- 
plification length the resonant character of the scattering 
is even more pronounced and £ osc of course also strongly 
depends on the gain value. The magnitude is even sig- 
nificantly larger than the amplification length £ a . This 
fact may represent a large obstacle in experiments. The 
physical picture presents itself now as the following, if one 
measures the intensity distribution between two points in 
coordinate space in the sample at a given distance r, i.e. 
at a finite value of Q in Eq. (|2"5|) . one will measure a dif- 
fusive intensity decay as characterized by the real part of 
the diffusion constant plus an exponential amplification 
due to £ a and additionally there is an amplitude modifi- 
cation proportional to cos{2-K£ a / £ osc ) in case one chooses 
the measuring distance to be £ a . This means, in general 
even though the intensity has experienced an exponential 
increase by a factor of e, this oscillating modulation fac- 
tor of the intensity is close to unity and therfore possibly 
hard to detect. 

Slightly changing the point of view let us now turn to 
the time scales. In the upper panel of Fig. [5] we present 
the gain induced growth time r a defined in Eq. (|40|) . 
For obvious reasons, r a displays a qualitatively similar 
behavior as the above discussed £ a , i.e., for frequencies 
within the scattering (Mie) resonances, the time to ex- 
ponentially increase the intensity is significantly smaller 
than for frequencies outside this range. Which reflects 
the fact that the gain coefficient is confined to the scat- 
terers volume only. 
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FIG. 7: Displayed in the upper the two particle amplification 
length l a for different values of the optical gain as a function 
of the dimensionless frequency. 

Displayed in the lower panel is the two particle oscillation 
length l for different values of the optical gain as a function 
of the dimensionless frequency. 



— e =10- le-4I 




FIG. 8: Displayed in the upper panel is the two particle am- 
plification time T a for different values of the optical gain as a 
function of the dimensionless frequency. 

In the lower panel we show the critical length Rcrit, character- 
izing the threshold volume in which the intensity experiences 
a laser-like growth behavior. 



Using the gain induced growth rate r a as defined in 
Eq. (|4T)]) , the intensity Green's function Eq. (f2"B1 ) may 
now be rewritten as 

a 

P<y ® , ™ ~~ -itt + iQ 2 lmD + Q 2 RcD - 1 /r Q ^ 

where the coefficient a may symbolically contain all the 
factors explicitly shown and discussed in Eq. (|28p. By 
inspection of the above equation, Eq. (|42j) and compar- 
ison with the Green's function Eq. (|35|) . it is to be rec- 
ognized that the energy density P(Q, f2) also exhibits a 
laser-like threshold behavior in complete analogy to the 
single-particle Green's function. 

Before discussing this threshold behavior in detail, we 
want to remind that our theory started with calculating 
the electrical field-field-correlator at different positions 
and frequencies Eq. ^ eventually leading to the evalua- 
tion the two particle Green's function given in Eq. (|28|) . 
This means, the momentum Q appearing in Eq. (|28j) rep- 
resents in Fourier space this relative position within the 
sample. In three dimensions the momentum Q therefore 
defines a volume unit within the sample. This volume is 
carefully to be distinguished from all over length scales 
e.g. the sample volume etc. It is merely the volume 
within one considers correlation effects of the diffusing 
behavior of the intensity. 

In analogy to ImE = in the single particle Green's 
function the threshold condition for the energy density 
now reads as follows 

Q 2 ReD - l/r a > (43) 
4tt 2 

^ ReJ D-l/ Ta = (44) 

**-crit 



leading to the definition of a third, a critical, length scale 
R crlt = 2ir^T a RcD. (45) 

This length describes the volume over which the energy 
density or intensity can compensate the diffusive loss by 
amplification due to the finite optical gain. 

As in the above described case of a single particle 
Green's function, the value R cr it marks the point in pa- 
rameter space where the laser threshold has been crossed. 
From the theoretically and experimentally known behav- 
ior of a lasing Mie sphere&22r— and the theoretical de- 
scription of the self-energy E(w) by the single scatterer 
t- matri x 32 ' 33 , and the definitions in Eq. (|40|) and Eq. 
(f3T)|) it becomes clear that in the limit of reaching the 
laser threshold within a Mie resonance r a is approaching 
zero, r a — > 0. And therefore the corresponding critical 
volume of the light intensity becomes point-like. In situ- 
ations where the optical gain is still below its threshold 
value with respect to the single scatterer Mie resonance, 
there is consequently a finite r a and therefore a finite crit- 
ical volume described by R cr it ■ In the limit of large gain 
it becomes clear from the definitions in Eq. ([40|) and Eq. 
([551 that R cr it in Eq. (|4"5"j) approaches the same order 
of magnitude as £ a - Which is perfectly meaningful, since 
in this limit the smallest and all-dominant length scale is 
set by the amplification length. For the most interesting 
intermediate range, we show the critical length R cr it in 
Fig. [5] as a function of light frequency for different gain 
strengths. Already for the gain values discussed in this 
publication, which are significantly below threshold, the 
ratio Rcrit I becomes as small as approximately xxx, 
c.f. lower panels of Figs. 181 andfTUl 

Finally, we emphasize that, once a Mie resonance is 
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FIG. 9: The imaginary part of the self-energy in the vicinity 
of the fifth Mie resonance for different but fixed values of gain. 
Increasing gain clearly yields narrowing and deepening until 
the laser threshold is reached, the curve approaches the shape 
of a delta-function. Beyond this point the theory is not valid 
anymore due to physical reasons. For completeness we also 
show the last curve Ime s = —0.0026 representing a system 
which has obviously crossed the threshold. The vertical line 
marks the frequency discussed in Fig. 1101 see also text below. 

close to lasing or the gain is very strong the growth time 
r a may become small and therefore the critical distance 
Rent may also become small as shown in the lower panel 
of Fig. [8j In general the length R cr it is not restricted 
to values above e.g. the single particle scattering mean 
free path. This is because the underlying physics is not 
scattering but frequency independent amplification. 

C. Gain and Laser Threshold 

In this last subsection, we want to closer discuss the ap- 
pearance of a laser threshold within our theory. Since we 
approximate the single particle self-energy by the single 
particle scattering matrix calculated within Mie theory, 
we first recall some basic and well known facts. The reso- 
nant features representing the resonant scattering modes, 
arise due to poles in the scattering coefficients, forming 
the t-matrix. These poles, or zeros of the coefficients' de- 
nominators, occur at complex frequencies, the closer the 
pole happens to be to the real frequency axes the more 
pronounced is the feature, i.e. the resonance becomes 
narrower and deeper. The effect of optical gain modeled 
as an imaginary part of the dielectric function, as it is 
done in this paper, is to lift the poles, i.e. it shifts the 
complex poles towards the real axes. In this way a gain 
narrowing is observed. The gain value corresponding to 
infinitesimal width of the resonance is believed to corre- 
spond to the experimentally observable laser threshold 8 . 
This situation coincides with a scattering pole right on 
the real frequency axes, i.e. the scattering resonance ap- 
proaches the limiting shape of Dirac's delta function with 



FIG. 10: The upper panel shows the imaginary part of 
the self-energy at a fixed frequency (cj/ojq — 2.1710269) as 
a function of gain, c.f. Fig. [9] The vertical line signal- 
izes the zero of ImS. From Fig. [9] and the text below, it 
is clear that this zero indicated that the laser threshold has 
already been crossed. The points (a) through (d) correspond 
to -Ime scat = {0.01,0.018,0.022,0.026}, i.e. to the curves 
shown in Fig. [9] The intersections with the vertical line there 
yield the values shown here. 

The lower panel displays the critical length i? crit as a function 
of gain. Due to the off-resonant frequency, R CTit has a local 
minimum and increases again with increasing gain. 



negative sign. Using gain values above threshold causes 
poles in the upper complex frequency half-plane, consti- 
tuting unphysical behavior due to the neglected laser dy- 
namics. If the complex frequency poles are in the upper 
half-plane, the theory still predicts a resonance feature, 
just with the "wrong" sign, i.e. instead of dips in the re- 
tarded self-energy one now observes peaks, see also Fig. 
[9] Additionally the larger the gain, the less pronounced 
the feature becomes, because the poles are then pushed 
away from the real axes in the complex frequency plane. 
This contains the risk of utilizing much to high gain val- 
ues, with poles in the negative complex frequency half- 
plane being at large distances to the real axes and there- 
fore yielding very weak features that do not necessarily 
break causality of the single particle Green's function for 
instance and might be overlooked at first glance. 

To illustrate the above discussed subject in and out of 
resonance, we considered a system slightly off-resonant 
but close to the fifth Mie resonance as depicted in Figs. 
and [TD] In Fig. [5] the imaginary part of the self en- 
ergy as a function of frequency is presented, the quantity 
describing both single particle scattering and amplifica- 
tion, as discussed earlier. For different gain values the 
narrowing and deepening of the resonance is clearly vis- 
ible as well as the difference of systems below and above 
threshold. The vertical (magenta) line marks the single, 
off-resonant frequency for which we study the self-energy 
(upper panel) and the critical length R cr it (lower panel) 
as a function of increasing gain, as shown in Fig. 1101 
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The behavior of Im£ is easily understood by compari- 
son with Fig. IH1 It is to be noted that neither the local 
minimum nor the zero define the laser threshold. The 
behavior of R cr it as defined in Eq. (|4"5)) is then shown 
in the lower panel of Fig. [10] The final increase is a 
consequence from the gain narrowing of the resonance. 
This decreases both the single particle scattering rate 
and the intensity emission rate l/r a defined in Eq. (|40[) 
and therefore increasing the critical volume. Following 
the above given line of arguments, it becomes clear that 
once a frequency closer to the resonance is chosen, the 
minimum value of R cr it decreases because the scattering 
rate and the intensity emission rate both increase. Even 
for a finite spectral width of the resonance, i.e. for a sys- 
tem below threshold, the value of the critical length may 
become quite small. 

IV. CONCLUSION 

In conclusion, we have presented a semi-analytical the- 
ory for scalar waves propagating in random, dissipating, 
i.e. also amplifying, media. The focus has been put on 
the influence of localization effects and finite gain on in- 



tensity transport in general. We found that for reason- 
able magnitudes of gain, the impact on transport quan- 
tities as e.g. the real part of the diffusion constant, scat- 
tering mean free path etc. is rather small. However, 
the gain introduces three new length scales natural to 
such systems, an amplification length £ a , an oscillation 
length £ osc and a critical length R cr it- The latter de- 
scribing the critical volume in which the intensity expe- 
riences a laser-like threshold behavior. The former two 
length scales constitute a growth length competing the 
diffusive loss of the intensity and the oscillation period of 
the intensity, respectively. Due to its comparably large 
magnitude, this oscillation length might be difficult to 
detect in experiments. We point out, that the critical 
length, or equivalently volume, has no lower bound other 
than zero since, e.g. for a Mie resonance reaching its las- 
ing threshold this volume becomes point like. For cases 
below threshold gain within the Mie scatterers or off- 
resonant light, the critical volume is finite and strongly 
influenced by the gain coefficient. 
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